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Abstract 

We attempt to generalize the familiar covariantly conserved Bel-Robinson tensor B^^ ~ 
RR of GR and its recent topologically massive third derivative order counterpart B ~ R DR, 
to quadratic curvature actions. Two very different models of current interest are examined: 
fourth order D = 3 "new massive" , and second order D > 4 Lanczos-Lovelock, gravity. On 
dimensional grounds, the candidates here become B ~ DRDR + RRR. For the D = 3 model, 
there indeed exist conserved B ~ dRdR in the linearized limit. However, despite a plethora 
of available cubic terms, B cannot be extended to the full theory. The D > 4 models are not 
even linearizable about flat space, since their field equations are quadratic in curvature; they 
also have no viable B, a fact that persists even if one includes cosmological or Einstein terms 
to allow linearization about the resulting dS vacua. These results are an unexpected, if hardly 
unique, example of linearization instability. 



1 Introduction 



The Bel-Robinson tensor B^ap was first discovered in ordinary Einstein gravity (GR) at D = 4, 
in a search for a gravitational counterpart of the usual matter stress-tensor T^ v . Since there can 
be neither local tensorial gravitational candidates of second derivative order (because quantities 
~ dg^,u dg a p, being frame-dependent, can be made to vanish at any point), nor any non-covariantly 
conserved ones, the successful candidate was indeed quadratic in the tensorial "field strengths" - 
curvatures B ~ RR, and covariantly conserved, in analogy with the quadratic Maxwell Tn V ~ F F. 
Subsequently, conserved B were found for arbitrary dimension and matter analogs have also been 
constructed (see, e.g., [1] for earlier references). 

Very recently, we showed that a conserved B ~ RDR could be defined for topologically 
massive gravity [2]. Given this theory's third derivative order, we speculated there on extending B 
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to other gravitational systems, in particular to the currently popular quadratic curvature models. 
The present note reports a negative outcome in two active, very different theories. For the fourth 
derivative D = 3 new massive gravity (NMG) [3], without an Einstein term [4] for simplicity, 
there is a B ~ dRdR at linearized level but it cannot be extended to the full theory, despite 
an enormous number of possible cubic correction terms. The Lanczos-Lovelock (LL) [5] models' 
quadratic curvature field equations, ~ RR = 0, obviously cannot even be linearized about flat 
space and turn out to have no conserved B either. Adding cosmological or Einstein terms to 
permit linearization (to effective cosmological GR) about the resulting dS vacua only allows the 
usual linear B ~ RR of GR, but no nonlinear extension. 



2 Fourth Order Models 



Since the machinery is considerably more complicated here than in GR, we analyze the (purely 
quadratic part of) NMG [3]: 

/ = J d 3 x ^ (R^ R^ - R 2 ) , = R^u-\ 9»v R- (1) 

[The Cotton tensor of TMG is the curl of the Schouten tensor R^.] The resulting field equations, 

DR^ - D^D V R — 4 R° R vc + RR^ V + R 2 = (2) 

of course obey Bianchi identities since (1) is an invariant. The field equations begin as DDR ~ 0, 
so B cannot just imitate the RR of GR, but must instead start as B ~ DRDR, along with possible 
cubic, RRR, terms since they have the same dimension as DR DR. Schematically then, we expect 
that 

B ~ DRDR + RRR. (3) 

[One can also consider similar terms ~ D(RDR) or RDDR, but they are already irrelevant at 
linear level.] We begin with the linearized (about flat space) truncation, ddR = 0, of (2); the 
corresponding linearized B is simply given by 

B^va/3 = Raf3° {R<rv,fM ~ Rfiv.a) ■ (4) 

Its conservation, 

0"£^ aj 8 = O (5) 

is verified using the field equations, the Bianchi identity, and antisymmetry of the parenthesis in 
(4). Conservation also holds for all permutations of i?'s (uaf3) indices; for completeness, we note 
also the identically conserved, if irrelevant, b^ap = D 1 H^^ va p, H antisymmetric in 7 ■<->■ fi. 

Thus encouraged, we turn to candidates (3) at nonlinear level, where the now covariant 
derivatives no longer commute. At first sight, there are so many available independent terms 
cubic in R^ v and g^ u R, that success seems guaranteed; however, it turns out to be unattainable. 
The simplest procedure is to take the covariant divergence of the initial B ~ DR (DR — DR) ansatz 
and try to compensate for the resulting cubic ~ RRDR terms by adding suitable RRR to B. The 
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procedure is straightforward, using the D = 3 properties 



Rpvpcr — 9 ftp Rvo Qvo Rpp Qvp Rp,a 9 per Ri/pi 

[D a , Dp] = \fi R aa - 51 R aP + g aa R\ - g a0 R\] V a . 



We find 



(6) 



D^B 



pvotfi 



R Pfl Rp/3 {Rpu;a — Rav;p) + ^ Ra/3;a {R UP Rpv ~ R° v R) 



+ (a 0) , (7) 



to which could be added the above (a/3is) permutations. The problem is now to find, among all 
(RRR)^ ua /3, those whose //-divergence cancels (7), or the sum of permutations. Note that (7) 
contains both fi- and a , /3-derivatives, and that the latter must arise from cubic terms of the form 
~ 9pa(R R R)v/3i an d possibly g^ a g u p(RRR); these are also plentiful. This systematic approach 
failed: any compensating term also created a new problematic one. We then undertook the brute 
force approach, involving the sum of all dimensionally possible candidate terms, ~ DRDR + 
R DDR + RRR; to be sure, any DR DR combinations that fail the linearized conservation test 
could be excluded a priori. 

We sketch the procedure and counts, but omit the uninstructive, laborious, detaila^- Start 
with the thirty nine possible terms of the form RRR, incuding g RRR terms, together with the 75 
D(RDR) candidates, giving a total of m = 114 potential terms. Take the divergence of each of 
these, enforcing the on-shell, Bianchi and three-dimensional simplifications. That gives, for each 
term, a right-hand side that is a linear combination of (DDR) DR, RRDR and R DDDR. There 
are n = 231 unique terms of this form in these sums. We can view the action of D^ as a linear 
operator taking the m starting ingredients to a space of n outputs: assign to each candidate term 
a basis vector {e^}™ 1 , and to each term appearing in the sums e» 3 the basis vector {Ej}" =1 , the 
matrix D £ K nxm j us t maps to its sum: Be^ = Y17=i ai We- constructed this matrix using a 
symbolic algebra package, and found that O has no null space, and so no conserved Bel-Robinson 
tensor exists. 



3 Second Order, LL, Theories 



At the other end of the quadratic curvature spectrum from the D = 3 model (1) lie the D > 4 LL 
theories, whose field equations remain of second derivative order, but do not admit any linearized 
expansion about flat space. The Lagrangians are, in vielbein formulation, 

r _ nbcdf... ABCDF ... r> R (o\ 

L — e e riabAB K-cdCD CfF • • • \°) 

The field equations simply consist of removing one vielbein factor, 

E fF = (eeRRe) fF = 0, (9) 
1 These may be found in the notebook at: http://people.reed.edu/~jfrankli/NullSpace.nb. 
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since the curvatures' variations vanish by the cyclic identities. Lowercase indices are local, capitals 
are world. Diffeomorphism invariance of the action ensures that DpE^ F vanishes identically, by 
the cyclic identities. Indeed, this is even true in D > 5: if we open any second pair of indices (gG), 
there results an identically conserved, hence uninteresting, fff9- FG ... We cannot even start with 
a linearized B, since the field equations (9) are intrinsically nonlinear. That there is no true B is 
rather clear: for example, attaching a curvature to H above loses conservation, which no additional 
cubics can help, while no DRDR term is relevant just because it is linearizable. Extending LL to 
include a cosmological or Einstein term does allow for linearization about a nonflat, dS, vacuum: 
a dS metric ansatz reduces (9) to the equation ~ a A 2 = 6 A in either choice. There is always 
a flat, A = 0, but also now a dS, A = b/a, vacuurnE). About the latter, the system linearizes to 
cosmological Einstein, which of course enjoys the same conserved B ~ RR as A = GR. However, 
just like its predecessor in NMG, this B evaporates as soon as we go to the full theory, since there 
is still no way to take into account the nonlinear RR part of (9). 

4 Conclusion 

We conclude that both NMG and Einstein+LL violate Fermi's rule: if the lowest order works, the 
conjecture is correct, while pure LL doesn't even have a lowest (here linearized) order. Clearly 
there is no invariance here like that [1] underlying ^-conservation in GR. 
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References 

[1] S. Deser, in "Gravitation and Relativity in General" (ed F. Atrio-Barandela and J. Martins), 
World Publishing (1999), |gr-qc/9901007l S. Deser and Z. Yang, Class. Quant. Grav. 7 1491 
(1990). 

[2] S. Deser and J. Franklin, Class. Quant. Grav. 28 032002 (2011). larXrv:1011.42401 

[3] E. A. Bergshoeff, O. Hohm and P. K. Townsend, Phys. Rev. Lett. 102 201301 (2009), hep- 
th/0901.1766. 

[4] S. Deser, Phys. Rev. Lett. 103 101302 (2009), hep-th/0904.4473. 

[5] C. Lanczos, Ann. Math. 39 842 (1938); D. Lovelock, J. Math. Phys. 12 498 (1971). 

[6] D.G. Boulware and S. Deser, Phys. Rev. Lett. 55 2656 (1985). 



While the appearance of a second, dS, vacuum in 72+LL - indeed in all R + R 2 models but (D = 4) i?+Weyl 2 - 
is reasonably clear [6], it may seem counter-intuitive that, unlike GR, quadratic models with cosmological term allow 
both flat and dS vacua. 
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